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Abstract 

Phase-space Lagrangian dynamics in ideal fluids (i.e, continua) is usually related to 
the so-called ideal tracer particles. The latter, which can in principle be permitted 
to have arbitrary initial velocities, are understood as particles of infinitesimal size 
which do not produce significant perturbations of the fluid and do not interact 
among themselves. An unsolved theoretical problem is the correct definition of their 
dynamics in ideal fluids. The issue is relevant in order to exhibit the connection 
between fluid dynamics and the classical dynamical system, underlying a prescribed 
fluid system, which uniquely generates its time-evolution. 

The goal of this paper is to show that the tracer-particle dynamics can be exactly 
established for an arbitrary incompressible fluid uniquely based on the construction 
of an inverse kinetic theory (IKT) (Tessarotto et al., 2000-2008). As an example, 
the case of an incompressible Newtonian thermofluid is here considered. 
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1 Introduction 



A basic aspect of fluid dynamics is related to the deflnition of the Lagrangian 
dynamics which characterizes both compressible and incompressible fluids. 
The customary approach to the Lagrangian formulation is based typically on a 
conflguration-space description, i.e., on the introduction of the (conflgur at ion- 
space) Lagrangian path, r(t), spanning the conflguration space (fluid domain) 
fi. Here r(t) denotes the solution of the initial- value problem ^ = V(r,t), 
with r{to) = Tq. Here = ^ -|- V ■ V denotes the so-called "fluid" convective 
derivative, Tq an arbitrary vector belonging to fl and V(r, t) the velocity fluid 
fleld, to be assumed continuous in f2 (closure of f2) and suitably smooth in f2. 
However, in turbulence theory the statistical formulation for the associated 
joint probability density for velocity increments requires the introduction of 
a phase-space representation of suitable type [H, in which usually the phase- 
space is identifled with F = i7 x R'^ (with closure F = i7 x R^), here denoted 
as restricted phase-space. Therefore it is natural to seek possible phase-space 
representations of this type for fluid systems. The goal of this investigation is 
concerned with the formulation of restricted-phase-space Lagrangian dynam- 
ics in such a way that the phase-space F coincides with the direct product 
space T = Q X V, Q being the fluid domain and V (velocity space) the set R'^. 
In particular in this paper, extending the formulationpreviously developed 
(Cremaschini et al. [2| and Marco Tessarotto et al. |3|, |j]), we shall adopt 
for this purpose a so-called phase-space inverse kinetic theory (IKT) (see also 
Tessarotto et al, 2000-2008 0, 0, 0, S [ij ) . Basic feature (of such an ap- 
proach) it that is relies on flrst principles, i.e., classical statistical mechanics 
and a prescribed complete set of fluid equation. This permits us to advance 
in time the relevant fluid flelds by means of phase-space Lagrangian equations 
deflned by the vector fleld X(x, t), namely 



where Xq is an arbitrary initial state of F (closure of the phase-space F). The re- 
sult appears relevant in particular for the following reasons: 1) the Lagrangian 
dynamics here determined permits to advance in time self-consistently the 
fluid flelds, i.e., in such a way that they satisfy identically the required set of 
fluid equations. For isothermal fluids, this conclusion is consistent with the re- 
sults indicated previously 0; 2) the Lagrangian dynamics takes into account 
the speciflc form of the phase-space distribution function which advances in 
time the fluid flelds; 3) the theory permits an exact description of the motion 
of those particles immersed in the fluid which follow the Lagrangian dynamics 
(classical molecules). 
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Phase-space Lagrangian dynamics and particle dynamics in ideal fluids are 
closely related issues. In fact both must be uniquely described via a suitable 
complete set of fluid flelds {Z{r, t)} which deflne the fluid state. This refers, in 
particular, to the so-called ideal tracer particles, for which both self-interaction 
produced by the perturbations of the fluid flelds generated by the same par- 
ticles and binary collisions among them are negligible (in this sense they can 
therefore be intended also as " coUisionless" ) . It is well known, however, that in 
customary approaches (see for example Maxey and Riley, 1982 the equa- 
tions of motion for ideal tracer particle are only known in some approximate 
sense and therefore do not reproduce exactly the correct fluid dynamics. 

The purpose of this paper is to show that an exact solution can be reached 
for phase-space Lagrangian dynamics, and in particular for the conventional 
tracer dynamics, based on the formulation of a suitable IKT. By deflnition, an 
IKT must provide the complete set of fluid equations describing the fluid, by 
means of velocity moments of an appropriate phase-space probability density 
function (pdf). We intend to show that such a theory can be uniquely deter- 
mined in the framework of classical statistical mechanics by invoking suitable 
statistical assumptions on the IKT. In particular, we present here a theory 
which applies to incompressible Newtonian fluids, including both isothermal 
and non-isothermal fluids 0, In the following we intend to show that cus- 
tomary tracer-dynamics equations due to several authors - including Tchen 
(1947 fll), Corrsin and Lumley (1956_0), Buevich (1966 [li]) and Riley 
(1971 [l5|) and Maxey and Riley (1982 [ll|) - are incompatible with the exact 
phase-space Lagrangian formulation here obtained. 

In detail the plan of the paper is as follows. First, in Section 2 previous ap- 
proaches to ideal tracer particle dynamics are summarized. Second in Sec. 3 a 
comparison between Eulerian and Lagrangian phase-space approaches is pro- 
vided. Furthermore, in Sec. 4 the IKT for incompressible thermofluids is pre- 
sented. This permits us to determine the appropriate form of the vector fleld 
X(x, t). Next, in Section 5 the Lagrangian formulation of IKT is discussed in 
detail. The new set of phase-space Lagrangian equations are shown to advance 
uniquely in time the relevant fluid flelds of an incompressible thermofluid. As 
a basic consequence, in Section 6 we will derive the exact dynamics of ideal 
tracer particles (see below for deflnition), comparing it with previous results. 



2 Previous approaches to tracer particle dynamics 

The motion of small particles (such as solid particles or droplets, commonly 
found in natural phenomena and industrial applications) which can be in- 
jected in a fluid with arbitrary initial velocity, in practice, may be very dif- 
ferent from that of the fluid. The accurate description of particle dynamics. 



3 



as they are pushed along erratic trajectories by binary colhsions (in real flu- 
ids) and by fluctuations of the fluid fields (in ideal fluids), is fundamental to 
transport and mixing in turbulence [1]. It is essential, for example, in com- 



bustion processes [16 



17 



as 



in the industrial production of nanoparticles 
well as in atmospheric transport, cloud formation and air-quality monitoring 

The Lagrangian approach - denoted as Lagrangian 



of the atmosphere [18|, 119 



turbulence (LT) - has been fruitful in advancing the understanding of the 
anomalous statistical properties of turbulent flows 20|. In particular, the dy- 
namics of particle trajectories has been used successfully to describe mixing 
and transport in turbulence (igI . 21|. Nevertheless, issues of fundamental im- 
portance remain unresolved (see for example Refs. [i^, [i^ for recent results 
regarding the Lagrangian view of passive scalar turbulence). In the past, the 
treatment of Lagrangian dynamics in turbulence was based on stochastic mod- 
els of various nature, pioneered by the meteorologist Richardson 2^ (see also 
22, 2^). These models, which are based on tools borrowed from the study of 



random dynamical systems, typically rely - however - on experimental verifica- 
tion rather than on first principles. However, in most cases there remains a lack 
of experimental data to verify the reliability of such models j22|. Verification 
can be based, in particular, on the measurement of fluid particle trajectories, 
obtained by seeding a turbulent flow with a small number of tracer particles 
and following their motions with an imaging system. On the other hand, 
the accurate evaluation of the Lagrangian velocity in laboratory turbulence 
experiments requires measurements of positions of tracer particle by using a 
suitable tracking system able to resolve very short time (and spatial) scales. 
In practice this can be a very challenging task since particle motions must be 
measured on very short time scales. 

As for the theory itself, rigorous results have been scanty, probably because 
of the subject complexity. In the case of ideal or dilute real fluids, however, 
particle motion is necessarily coUisionless (in the sense specified above) while 
the dynamics of ideal tracer particles is controlled by the force produced on 
them only by the unperturbed fluid fields. Several authors have tried in the 
past to derive, based on phenomenological arguments, an approximate equa- 
tion for the ideal tracer-particle dynamics, describing the motion of a particle 
suspended in a non-uniform flow. Since the original Basset-Boussinesq-Oseen 
(BBO) equation 25|, l26|, |27|], formulated in the case of a uniform flow, sev- 
eral papers have appeared proposing modifications or corrections of the same 
equation for non- uniform flows (for a review see [llQ. The first attempt at a 
generalization of this type is due to Tchen, (1947 12]), who considered the 
motion of a rigid sphere in an incompressible isothermal Navier-Stokes (NS) 
fluid. Tchen 12j derived an approximate equation of motion for a finite-size 
spherical particle of radius a and mass mp describing its dynamics in terms 
of the Newtonian state of its center of mass {r{t),\{t)}. His equation was 
later modified by Corrsin and Lumley (1956 ^]), to take into account con- 
tributions due to pressure gradients previously ignored, and by Buevich (1966 
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14l |) in order to consider also the effect of viscous stress. The version of the 
equation currently adopted by some authors (see for example Gui et ai, 2008 
"iq where it was used to investigate modifications of turbulence) is, however, 



the one later developed by Maxey and Riley (1982 [ll|) in which also the 
buoyancy contribution produced by the volume displaced by the particle was 
taken into account. In the approximation in which perturbations of the fluid 
fields produced by the particle are negligible the equation of motion developed 
by Maxey and Riley reduces to: 



d L)V(x,t) 
mp—v{t) = mp 



dt ^ ' Dt 
{mp - mp) g 



x=r(t) 



^"^F^{v(t)-V(x(t),t)} 



f2) 

x=r(t) 



{M-R equation), where the last contribution on the r.h.s. denotes the so-called 
buoyancy effect. Here the notation is standard. Thus, {V(r, t),p(r, t)} are re- 
spectively the fluid velocity and pressure, f the volume force density and finally 
Poji' > the constant mass density and kinematic viscosity (with u related 
to the dynamic viscosity by the identity u =/i/po)- In particular, f can be 
written f = — V0 -|- f^, with pog= — V0, where g is the constant local gravity 
acceleration, (p = pogz the gravitational potential (hydrostatic pressure) and 
fp a possible non-potential force density. Moreover, mp is the mass of the fluid 
displaced by the sphere, ^ = ^ + v(t) ■ V the "particle" convective deriva- 
tive and V(r(t),t) and v(t) — V(r(t),t) are respectively the fluid and particle 
relative velocities evaluated at the position of the particle center of mass. It 
should be noted, however, that also this equation is still unsatisfactory. In 
fact, it is obtained by requiring that the particle velocity remains always suit- 
ably close to the fluid one, so that contributions due to the relative velocity - 
i.e., proportional to the particle relative velocity u(t) = v(t) — V(r(t), t) - are 
actually ignored in Eq.(l2]), requiring that for all {r(t),v(t)} there results: 

|u(t)|«|V(r(t),t)|. (3) 

The limitation appears serious because the actual dynamics (i.e., both the 
velocity and the acceleration) of ideal tracer particles may be in principle very 
different from that of the fluid elements. In addition, the accurate descrip- 
tion of ideal particle dynamics is essential both in LR and in environmental 
fluid dynamics (dynamics of anthropogenic pollutants in the atmosphere, dif- 
fusion of dusty particles, droplets, aerosol particles, etc.). This involves, in 
fact, the ability to simulate tracer dynamics in a variety of different physical 
conditions and in fluid flows characterized by a turbulent behavior. Therefore, 
an open issue remains the very definition of the dynamics of tracer parti- 
cles which may be injected in an ideal fluid with arbitrary initial velocities. 
Clearly, such a formulation - if achievable at all - should rely exclusively on 
first principles, i.e., in particular the exact validity of the fluid equations. This 
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problem is closely related to the formulation of phase-space approaches for 
fluid systems, based on the introduction of suitable phase-space representa- 
tions of classical fluid dynamics in terms of an appropriate phase-space pdf. 
By construction, the dynamics of a fluid is completely described by the time- 
evolution of its fluid fields (which, in turn, are assumed as classical solutions 
of a well-posed initial-boundary value problem defined by a complete set of 
fluid equations). Therefore, unless suitable restrictions are posed, phase-space 
dynamics remains intrinsically non-unique. This is true, in particular, due to 
the arbitrariness in the choice of the possible phase-space and the definition 
of the evolution equation for the pdf. 



3 Eulerian and Lagrangian phase-space approaches 

As it is well known, phase-space descriptions of fluids can be achieved in 
principle choosing either an Eulerian or a Lagrangian point of view. Based on 
the IKT approach for incompressible fluids earlier developed j^, 0, Q, @, 0] such 
a connection can be uniquely established. IKT is based on the identification 
of the complete set of fluid fields (which describe the fluid) with velocity- 
moments of a suitably-defined kinetic distribution function /(x,t). The pdf is 
assumed to satisfy the basic principles of classical statistical mechanics, which 
include in particular: 

(1) the principle of conservation of probability; 

(2) a suitable correspondence principle (Ellero et al., 2005 [3]) whereby ap- 
propriate (velocity-) moments of the pdf can be identified with the rele- 
vant fluid fields; 

(3) the principle of entropy maximization (PEM, Jaynes, 1957 |29l] ) : 

(4) and an entropic principle assurinqthat the statistical entropy cannot de- 
crease in time (Tessarotto, 2008 |3|). 

The first axiom implies that the pdf must satisfy a Liouville equation, i.e., 
a Vlasov-type inverse kinetic equation (IKE, Ellero et al. [3]). As it is well- 
known, this type of kinetic equation is, in fact, appropriate for the statistical 
description of particles subject solely to mean-field interactions. This is con- 
sistent with the assumption of an ideal fluid, i.e., a continuum in which the 
fluid elements are, by definition, subject only to mean-field interactions. In 
such the time-evolution of the pdf is determined by the Eulerian IKE 

Lf{^,t)=0. (4) 

Here /(x, t) denotes the Eulerian representation of the pdf, L is the streaming 
operator Lf = ^/ + ^ ■ {X(x, t)f} , X(x, t) = {v, F(x, t)} a suitably smooth 
vector field, while v and F(x, t) denote respectively the velocity and an appro- 
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priate " mean- field" acceleration vector field. The implications of the principle 
of maximum entropy 29|] and of the entropic principle, both involving the 
assumption that the Boltzmann- Shannon (B-S) entropy functional 



S{t) ^ 5(/(x,t)) = - I cix/(x,t) In /(x,t) (5) 

r 

exists, have been discussed elsewhere ji], [sf (see also Tessarotto et al, 2007 
[ol). In particular this provides a well-defined initial condition for the pdf, 
/(x,to) = fM{^,to) [see below Eq. flT^ ] and also a (generally non-unique) rep- 
resentation for the streaming operator L (Tessarotto et al, 2006 |8|). As a 
main consequence the same approach can in principle be used to determine in 
a rigorous way the Lagrangian formulation for arbitrary complex fluids. Al- 
though the choice of the phase-space F is in principle arbitrary, in the case of 
incompressible isothermal fluids, it is found |6| that it can always be reduced 
to the direct-product space T = Q xV {restricted phase-space) , where fl,V <^ 
M.^, Q is an open set denoted as configuration space of the fluid (fluid domain) 
and V is the velocity space. 



4 The case of an incompressible thermofiuid 



Let us consider for definiteness an incompressible, viscous and generally non- 
isentropic thermofiuid (which comprises as a particular case alsc) the treatment 
of incompressible isothermal fluids earlier developed in Ref. j7|). This is de- 
scribed by the fluid fields {Z} = {p > 0, V,p > 0, T > 0, St} , to be identified 
respectively with the mass density, the fluid velocity, pressure, temperature 
and thermodynamic entropy. In the open set Q they are assumed to satisfy 
the so-called incompressible Navier-Stokes- Fourier equations (INSFE), i.e.. 



P = po 


> 0, 




V- V 


= 0, 






Po 


[Vp 


Dt 


= X V'T + 






>o. 





where po is a constant, 

1 



(6) 
(7) 



+ ^ 1 + = K, (9) 

PoCp 



(10) 



Fh = -— [Vp - f] + z/V^V (11) 
Po 
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is the fluid acceleration and jji = ^ + V-V the convective derivative. These 
equations are assumed to satisfy a suitable initial-boundary value problem 
(INSFE problem) so that a smooth (strong) solution exists for the fluid fields 
{Z} . Here the notation is standard ^]. Thus, Eqs.([n])-([HD together are denoted 
incompressible Navier-Stokes equations (INSE) [with Eqs.([6]) and ([7j) repre- 
senting respectively the so-called incompressibility and isochoricity conditions 
and Eq. (IHl) the forced Navier-Stokes equation forced Navier-Stokes equation 
written in the Boussinesq approximation], while Eq.(l9]) is the Fourier equation. 
Finally ([T]) denotes the entropy inequality which defines the second principle 
of thermodynamics . As a consequence, in such a case the force density f reads 
f = Pog (1 ~ kp^T) + fi, where the first term represents the (temperature- 
dependent) gravitational force density, while the second one (fi) the action of 
a possible non-gravitational externally-produced force. Hence f can be written 
also as f = — V0 + fR,where = pogz and f/j = —pogkp^^T + fi denote respec- 
tively the gravitational potential (hydrostatic pressure) and the non-potential 
force density. Moreover, in Eq. (Q K and J are the quantities of heat generated 
per unit volume and unit time by all sources and, respectively, only by the 
external sources. In particular, the inequality (ITU]) defines the so-called 2nd 
principle for the thermodynamic entropy St- For its validity in the following 
we shall assume that there results everywhere mVLxI 



which defines a so-called externally heated thermofluid. In these equations 
g,kpQ,h',x and Cp are all real constants which denote respectively the local 
acceleration of gravity, the density thermal-dilatation coefficient, the kinematic 
viscosity, the thermometric conductivity and the specific heat at constant 
pressure. Thus, by taking the divergence of the N-S equation ([8]), there it 
follows the Poisson equation for the fluid pressure p, namely V^p = — poV ■ 
(V ■ VV) + V ■ f , with p to be assumed non negative and bounded in Q x I. 

4-1 Functional uniqueness of IKT 

Let us now assume that /(x,t) is a solution of the Eulerian kinetic equation (jll) 
defined in a suitable extended phase-space Fx/, where J C M is a suitable time 
interval. In such a case, we intend to show that /(x,/!:) [to be assume strictly 
positive] can be defined in such a way that the fluid fields V,pi and St can 
be identified with its velocity moments / (ivG(x,t)/(x,t), where respectively 



(^(x,^) = v,poM^/3, — Z?T,/(x,t), u = V — V(r,t) is the relative velocity and pi 
the kinetic pressure defined as: 




(12) 



(13) 
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Here po{t) (to be denotes as pseudo-pressure 



is an arbitrary strictly posi- 



tive and suitably smooth function defined in /. Moreover, mp > is a constant 
mass, whose value remains in principle arbitrary. In particular it can be iden- 
tified with the average mass of the molecules forming the fluid. Finally, the 
thermodynamic entropy St can be identified with the Shannon statistical en- 
tropy functional 5'(/(x,t)), provided the function po{t) is a suitably prescribed 
function and /(x,t) is strictly positive in the whole set F x /. To reach the 
proof, let us first show that, by suitable definition of the vector field F(x, t), a 
particular solution of the IKE (jl]) is delivered by the Maxwellian distribution 
function: 



where vth = Y2pi(r, t)/po is the thermal velocity driven by the kinetic pres- 
sure pi{r,t). Based on the results earlier obtained for ideal isothermal and 
incompressible fluids , [qI and incompressible thermofiuids [sl the follow- 
ing theorem is reached: 

Theorem - IKT formulation for INSFE 

Let us assume that the INSFE problem admits a smooth strong solution in 
Fx/, such that the inequality [W^) is fulfilled and the fluid fields {Z} belong 
to the "minimal functional setting" (see Ref. ^/j. Moreover, let us assume 
that in Q X I : 

1) the pdf f{'^,t) is suitably smooth, strictly positive and admits the velocity 
moments G(x,t) = 1, v, pof ^, PoU^; PoUu; thus, we denote in particular Q = 
PoJ d^vu^f and 11 = PoJ d^vuuf; 

2) the B-S entropy integral ^ exists in the time interval / CM; 

3) the inequality [W^) is assumed to hold; 

4) there results identically (correspondence principle): 




(14) 




(15) 




(16) 




(17) 



Srit) = 5(/(x,t)). 



(18) 
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Then it follows that: 

Ti) the local Maxwellian distribution fM^'^.t), defined by Eg. PTP, is a partic- 
ular solution of the inverse kinetic equation 



T2) the mean-field acceleration vector field F reads 

F(x,t;/)=Fo + Fi. 



(19) 



The functional form of the vector fields Fq, F;^ is determined uniquely by re- 
quiring that they depend only on the velocity moments indicated above. They 
read respectively: 



Fo(x,t;/) = - V-n-Vpi + f^j +D(x,t)+z/V2V, 
Po L - J 



(20) 



1 1 1 . 1 



Fi(x,t;/) = -u -A+-V-Q 2 

2 ipi pi Pi 



1 r. 



where 



v-n 




d 



A = ^^{po+p)-V- 



■Q + 



(21) 



D(x,t)=-{VV-u + u- VV} 



Dt Po 



f22) 



PoK _ D 

— = TTtPi; (23) 
mp Dt 



T3) for an arbitrary pdf /(x,t) fulfilling assumptions 1-3 equations / l73]) -/ f^) 
are fulfilled identically in Q x I. 

PROOF 



Let us, first, prove proposition Ti. For this purpose, let us assume that a 
strong solution of the INSFE problem exists which in the set Q x I satisfies 
identically Eqs.(l6])-(l9]). In such a case it is immediate to prove that fM{'^,t) 
is a particular solution of the inverse kinetic equation (jl]). This can be proved 
either: a) by direct substitution of / = /a/(x, t) in Eq.(jl]) (Proposition A); 
b) by direct evaluation of the velocity moments of the same equation for 
G{x,t) = 1, v,M^/3 (Proposition B). Regarding Proposition B, we notice that 
the first two moment equations coincide respectively with the isochoricity and 
Navier-Stokes equations [Eqs. ([6]) and ([8])]. Therefore, the third moment equa- 
tion delivers the Fourier equation [Eq.([9])]. The same proof (Proposition B) is 
straightforward also if / 7^ /m(x, t). This is reached again imposing the same 
constraint equation (ITSll on first velocity-moment of the distribution function 
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/. The proof of T2 and follows in the same way by direct evaluation of the 
velocity moments of the same equation for G(x,t) = 1, v, po|^, poU^, PoUU. 
Q.E.D. 



5 Lagrangian formulation of IKT 



The previous results permit us to formulate in a straightforward way also the 
equivalent Lagrangian form of IKE [see Eq.(jn])]. The Lagrangian formulation 
is achieved in two steps: a) by identifying a suitable dynamical system, which 
determines uniquely the time-evolution of the kinetic probability density pre- 
scribed by IKT. Its flow defines a family of phase-space trajectories, here 
denoted as phase-space Lagrangian paths (phase-space LP's); b) by proper 
parametrization in terms of these curves of the pdf and the inverse kinetic 
equation, the explicit solution of the initial-value problem defined by the in- 
verse kinetic equation (jlj) is determined. First, we notice that - in view of the 
previous theorem it is obvious that the phase-space LP's must be identified 
with the phase-space trajectories x(t) of a classical dynamical system 

Xo x(t) = Tt^to^o = x(xo, to, t) (24) 

(here denoted as INSFE dynamical system), with Tt^to the corresponding evolu- 
tion operator generated by the vector field X(x, t), to be prescribed according 
to the previous theorem. Hence, the initial-value problem ([T]) is realized by 
the equations 

d 



d 



|r(t)=v(t), 
v(t) = F(r(t),t;/), 



V . ^25) 

v{to) = To, 
v{to) = Vo, 

where the vector field F(r(t), t; /) is defined by Eq. f|T9|) . Here, by construction 
v(t) and F(r(t), t; /) are respectively the Lagrangian velocity and acceleration, 
both spanning the vector space M^. In particular, F(r(t), t; /), which is defined 
by Eqs. fl2(I]) - fl^ . and depends functionally on the kinetic probability density 
/(x, t), is the Lagrangian acceleration which corresponds to an arbitrary kinetic 
probability density /(x, t). From the theorem it follows that in the Lagrangian 
representation the kinetic equation {Lagrangian IKE) can be written in the 
form 



J(x(t), t)/(x(t), t) = /(x„, to) = foi^o) (26) 
where /o(xo) is a suitably smooth initial pdf and J(x(t),t) is the Jacobian 

dx{ 

ax 



J{x(t),t) = of the map Xq x(i(:) which is generated by Eq. (!25|) . It 
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follows that the Lagrangian equation (126|) ] is uniquely specified by the proper 
definition of a suitable family of phase-space LP's. Eq. (!26|) also provides the 
connection between Lagrangian and Eulerian viewpoints. In fact the Eulerian 
pdf, /(x, t), is simply obtained from Eq. fl26p by letting x = x(t) in the same 
equation. As a result, the Eulerian and Lagrangian formulations of IKT, and 
hence of the underlying moment (i.e., fiuid) equations, are manifestly equiva- 
lent. 



6 The exact dynamics of ideal tracer particles 



Let us now analyze in detail the equations of motion for ideal tracer particles 
immersed in an incompressible thermofiuid described by INSFE [Eqs.fl6|)- f|T0|) ] 
. First, it must be remarked that two types of forces can in principle be present: 
a) a volume force, acted by the fiuid (the same one which is responsible of the 
phase-space Lagrangian dynamics); b) particle- localized forces, such as the 
gravitational pull, acting directly on the tracer particle. In particular, regard- 
ing the first one, its specific form depends on the assumed pdf to be associated 
to the fiuid. As discussed elsewhere [sOi] its choice depends closely on the type 
of fiuid to be considered, i.e., deterministic or stochastic (as appropriate to de- 
scribe turbulent fiows). In particular, the position / = /jvf(x, t), with /^(x, t) 
defined by Eq. ffMl) is suitable for the description of deterministic, i.e., non- 
turbulent fiows. In such a case the appropriate form of the equation is obtained 
from Eqs. f|T9l) - fl23l) . Instead, the general case, in which one allows / 7^ /Af (x, t), 
is provided by Eqs. (dH) with (I2UD and mh. 

Let us now assume that rrip mp, mp denoting the mass of the displaced 
fiuid. In view of the IKT approach the equations of motion depend necessarily 
on the form of the pdf /(x, t), and hence describe in this sense the conditional 
phase-space dynamics. To construct the equation of motion for an ideal tracer 
particle of arbitrary mass, let us now assume, for definiteness, that the sole 
particle-localized force acting directly on the tracer particle is produced by 
the gravitational pull. The equation of motion for an ideal tracer particle of 
mass mp reads simply in such a case: 



nip 



d , , 



mr [F(x,«; /) - g] . (27) 



which describes the conditional dynamics of an ideal tracer particle immersed 
in an incompressible thermofiuid described by a pdf /(x, t). We stress that 
the form of the pdf depends on the specific assumptions made on the fiuid 
'3o| . The interpretation of this equation is as follows. The terms on the l.h.s. 



represent the "inertial" and gravitational forces acting on the tracer particle. 
Instead, all the terms on the r.h.s. represent the volume force acting responsible 
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for the phase-space Lagrangian motion. The physical interpretation of the 
various contributions appearing in the volume force [see Eqs. fll9l). (!20|) and 
fl21l) ] is made transparent by representing them in terms of the vector fields 
mirFo(x,t;/) and mpFi{-x,t; f), to be interpreted as mean- field forces. One 
obtains in fact, in particular: 



Fo(x,t;/)-g-^Vpi 



2 



It is immediate to prove that the terms in the first equation take into account 
the fluid and convective forces mpFu and mpY){^,t), with Yh and D(x,t) de- 
noting respectively the fluid acceleration fill I) and the convective term defined 
above [see Eq. fl22p ]. In a similar way one can show that, in case mp ^ mp, 
Eq. fl271) recovers also the buoyancy force denbsity {mp — mp) g pointed out 



by Maxey and Riley [11[. Finally, for the sake of comparison, let us consider 
the case of an isothermal fluid and require - consistent with Eq.Q - that lo- 
cally in the extended phase-space T x I condition holds at time t = to- 
The validity of previous tracer-particle equations [i.e., in particular the MR 
(Maxey-Riley) equation ([2])] requires that the inequality ^) holds for all times 
(t G /). In particular, one can prove that in such a limit our Eq. fl27l) agrees 
with the MR equation ([2]). In fact, it yields 



d DV(-K,t) 



dt ^ ' Dt 



+ mpFi - {mp - mp) g, (29) 

x=r(t) 



where, in validity of ([3]), it follows that the pressure mean- field force can be 
approximated as mpFi = —\mp ^ {v(t) — V(x(t),t)} . On the other 

hand, a serious objection to all of the previous tracer-dynamics equations 
is provided by the possible violation of the asymptotic condition which 
may not be uniformly fulfilled in the set I. This occurs, manifestly, if an initial 
condition of the type |u(to)| ~ \Y{T{to),to) \ is imposed on a tracer particle (this 
requirement is not physically unreasonable since, in principle, tracer particles 
might be injected in a fluid with arbitrary initial velocities). However, even if 
initially (at t = to) one requires the validity of ([3]) in general it may well be 
also that |u(t)| ~ |V(r(t),t)| at some later time {t > to). This can be achieved, 
for example, even imposing the initial condition |u(to)| = 0. The result is a 
consequence of Eq. fl7r|) . Indeed one can prove that, even imposing the initial 
condition |u(to)| = 0, generally |u(t)| ^ 0, with |u(t)| not satisfying In 
other words, ideal tracer particles having initially the same local velocity of the 
fluid may develop in time a finite relative velocity. This implies that, generally, 
the full exact tracer-dynamics equation, i.e., Eq. (!271) . should be used, instead 
of the asymptotic approximation indicated above [see Eq. (12^ ]. 
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7 Concluding remarks 



In this paper tlie pliase-space Lagrangian dynamics lias been determined as 
appropriate for an incompressible thermofiuid described by a suitable set of 
fluid equations [INSFE, see Eqs.(!6l)-( fT0l) ]. We have shown that, based on the 
formulation of a restricted phase-space inverse kinetic theory, the phase-space 
Lagrangian dynamics can be uniquely established. The governing equations 
which determine the phase-space Lagrangian trajectories (LP's) are found 
to depend functionally on the pdf [/(x, t)], to be uniquely associated to the 
fluid by means of the IKT here adopted. In particular, the theory permits to 
advance uniquely in time /(x, t) and in terms of the same pdf also the complete 
set of fluid fields which describe the fluid. This feature is of fundamental 
importance in turbulence theory (see (3o|). 

As a further consequence, the dynamics of ideal tracer particles (i.e., for which 
the perturbations of the fluid fields produced by the same particles are negligi- 
ble) is established. Remarkably, this result overcomes limitations of customary 
ideal tracer-dynamics equations [see in particular the Maxey and Riley equa- 
tion [11] given by Eq.([2])]. All of these equations are, actually, in disagreement 
with the present theory for finite particle relative velocities. The basic new 
result is represented by Eq. fl271) which describes, for arbitrary initial veloc- 
ity, the conditional dynamics of an ideal tracer particle in an incompressible 
thermofiuid described by a suitable pdf /. 
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